Many-site coherence revivals in the extended Bose-Hubbard model 
and the Gutzwiller approximation 
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We investigate the collapse and revival of first-order coherence in deep optical lattices when long- 
range interactions are turned on, and find that the first few revival peaks are strongly attenuated 
already for moderate values of the nearest-neighbor interaction coupling. It is shown that the 
conventionally employed Gutzwiller wavefunction, with only onsite-number dependence of the vari- 
ational amplitudes, leads to incorrect predictions for the collapse and revival oscillations within the 
extended Bose-Hubbard model. We provide a modified variant of the Gutzwiller ansatz, reproducing 
the analytically calculated time dependence of first-order coherence in the limit of zero tunneling. 
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I. INTRODUCTION 

The Hubbard model describes interacting particles in 
sufficiently deep periodic lattice potentials, so that the 
tight-binding approach applies, and the Wannier func- 
tions form the basis of the bosonic or fermionic field op- 
erator expansion [ll-Q ■ It was realized about a decade ago 
that in optical lattices, generated by counterpropagating 
laser beams, a clean and highly controllable implementa- 
tion of the Hubbard model can be obtained [J] . Initially, 
attempts to achieve a precise understanding of the col- 
lective behavior of atoms or molecules trapped in the lat- 
tice was concentrated on the classic version of the model 
with onsite interactions only. However, a zoo of quantum 
phases beyond the superfluid and Mott phases, like var- 
ious supersolid and charge-density wave phases, emerges 
when long-range interactions are turned on [5| . The most 
prominent and (for electrically neutral atoms/molecules) 
physically realizable examples of the latter are magnetic 
and electric dipole-dipole interactions [^]. The so-called 
extended Hubbard model which incorporates the long- 
range character of the interactions Q- It is crucial for 
the determination of possible phases to obtain informa- 
tion on the interaction parameter values in the extended 
Bose-Hubbard model, and we suggest, in the following, 
experimental means to measure them by using the tem- 
poral variations of first-order coherence. 

Collapse and revival oscillations are familiar from 
quantum optics [8| as a sensitive measure for the spec- 
trum of the Hamiltonian and the many-body dynam- 
ics of the system. In ultracold bosonic quantum gases, 
collapse and revival oscillations were experimentally ob- 
served both on the two- and multi-body interaction level 
[9|, |lO[, when switching quasi-instantaneously from the 
superfluid to the Mott side in a dynamical quantum 
phase transition, a quench [ll|. Rapid experimental 
progress in the meantime even allows for the single-site 
and single-atom resolution of the time evolution of the 
site occupation number [12|, ^^ . 

In the optical lattice realization of collapse and revival 
of coherence, effectively many small samples on the var- 



ious sites are prepared, which are independent of each 
other when only contact interactions are present. When 
the particle numbers in individual sites are small, i.e. at 
small filling, the revival times are short enou gh f or ob- 
serving the collapse and revival phenomenon [IJ]. The 
collapse and revival times have in addition been demon- 
strated to be sensitive to the effective dimension of a 
single well and its trap power law [l5|. 

It has been previously shown that a sensitive measure 
of the many-body dynamics due to collapse and revival 
dynamics in optical lattices may be experimentally ac- 
cessible when tunneling across different sites is included 
[l6H18t . Furthermore, significant effects on collapse and 
revival can be expected when effective three-body inter- 
actions [l9[ as well as the number squeezing of a initial 
state with finite interactions on the superfluid side are 
taken into account [20[. In addition, the number squeez- 
ing of an initial state can be caused by flnite sweep rates, 
when quenching from the superfluid to the localized Mott 
phase lU. 

We explore in the following effects which are intrinsi- 
cally beyond the familiar purely onsite collapse and re- 
vival evolution in optical lattices and come from long- 
range interactions (that is not those, e.g., coming from an 
additional external conflnement potential [23). The ex- 
tended Hubbard model provides nearest-neighbor (NN) 
coupling of sites on the two-body interaction level, so 
that the Hamiltonian converts the onsite phenomenon 
collapse and revival obtained for contact interactions into 
a genuine many-body-many-site effect. We demonstrate 
that an attenuation of the first few phase coherence re- 
vivals takes place in a range of the ratio of onsite and 
NN interaction couplings relevant to the predicted quan- 
tum phases of the extended Bose-Hubbard model [^, [2j1 ■ 
It is shown that, as opposed to the reduction of the co- 
herence signal by the (phase-sensitive) tunneling term in 
the Hamiltonian, no real damping takes place due to a 
NN interaction coupling, instead there are two distinct 
revival scales related to onsite and NN interaction cou- 
plings U and V, respectively. The consequent effective 
attenuation of the first few revival peaks due to V is 



much stronger than those to tunnehng damping, increas- 
ing with fining rather than decreasing, and can easily 
overcome inhomogenity (external trapping) efi^ects. Wc 
therefore propose a sensitive measure of the value of V 
realizable in current optical lattice experiments. 

The modulation of the first-order coherence signal with 
frequency V/h, taking as in Sec. II below a coherent 
product state as the initial state after a quench into a 
deep lattice, does not occur in the conventional time- 
dependent Gutzwiller approach. This approach employs 
a product state of purely onsite-oocupation dependent 
variational wavefunctions, and is indeed widely used for 
ground-state calculations in the extended Bose-Hubbard 
model, cf. e.g. (a. l24| - [2a | . As we will show below, it how- 
ever leads at least in dynamical situations (i.e. with ex- 
plicit time dependence of correlation functions) to in gen- 
eral grossly incorrect predictions and hence needs modifi- 
cation. We will devise a variant of the Gutzwiller product 
state with the variational amplitudes depending on the 
NN particle-numbers. The form of this wavefunction is 
based on the analytical evolution of first-order phase co- 
herence in deep optical lattices, i.e., when tunneling is 
exponentially small and thus negligible. We then apply 
this variant of the Gutzwiller wavefunction, appropriate 
to the extended Bose-Hubbard model, to an initial state 
which is number-squeezed. We again compare the re- 
sult with that obtained by the conventional Gutzwiller 
approach, showing that in most experimentally relevant 
cases the NN-coupling-induced attenuation of the coher- 
ence signal will completely dominate the damping caused 
by residual tunneling. 



II. COHERENT INITIAL STATE 



The Hamiltonian of the single-band extended Bose- 
Hubbard model is assumed to be of the form 



M M M 



where < ij > indicates summation over NN, and M is 
the number of lattice sites. We consider here the most 
general situation that the parameters J, U, V can be con- 
trolled independently of each other, e.g. by optical lattice 
depth and the ratio of contact and long-range interac- 
tions. Their values are determined by overlap integrals 
of Wannier functions incorporating single-particle energy 
operators and two-body interaction power laws and effec- 
tive ranges, respectively. The two-body interaction terms 
correspond to density-density type interactions with on- 
site ([/) and NN {V) parts. 

We take in this section the initial state to be the prod- 
uct of onsite coherent (superfluid) states with Poissonian 
probability distribution of site occupation 
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where \ni) are local number eigenstates, iii \ni) 
and di \ai) = at \ai) with average fii = {hi) = \a. 
a state is for example physically realizable with a quasi- 
instantaneous sweep across the superfluid-Mott transi- 
tion starting from (deep inside) the superfluid phase, in 
which J > [/,F @, [ll]. We note that finite M cor- 
rections due to a particle-number-conscrving initial state 
and the resulting time evolution were investigated for the 
Bose-Hubbard model in [18]. However, they became neg- 
ligible for a number of sites M > 5, so that we can expect 
to be able to safely neglect them for sufficiently large M 
also with nonvanishing off-site interaction coupling. 

The limits of either U and J being dominant over the 
other couplings arc well known, and lead to the (deep) 
Mott phase (at commensurate filling in the present ho- 
mogeneous case N/Ad = n integer, where TV is the total 
number of particles) and the superfiuid phase of weakly 
nonintcracting bosons, respectively. We choose [/ = 1 as 
the unit of energy and 27r/C/ as the unit of time [h = 1). 

We first derive an analytical result for the exact time 
evolution of first-order coherence when J = in ([T|) . Be- 
cause the energy is expressible in the number (i.e. Fock 
space) basis when J = 0, the first order correlation func- 
tion can be analytically determined. Calculating with 
the initial coherent product state ^ the off-diagonal el- 
ements of the single-particle density matrix, one obtains 
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where the site-occupation-dependent interaction energy reads £^{ni} — ^ X^i '^j("-i " -'^) + T X]<i7> "-j^jj *5,nd [il 



and [jl] indicate sums over all NNs / of i and j at fixed 
i and j, respectively. We assume for simplicity in ([3]) 
that |j — j| > 2, for which the coherence signal becomes 
independent of distance \i — j\ (expressions when \i — j\ < 
2 can easily be derived as well) . The last line is valid for 
homogeneous filling n. Finally, ly is the number of NN, 
^■^■Sfifci ~ ^' ^^'^ ^ij ~ 3''"g['^i] ~ 8'i'g[Q^i] describes a 
possible initial phase offset between sites. 






< 



1.0 




A 


A A « 


- 


0.8 








: 


0.6 






1 i1 1 


; 


0.4 






ft \\ 


A 


0.2 


1 


/ \ 


l/UAi 


\ 


00 


V. 


/ V^ 


' \J xJxJ 


/4=: 




3 


1 


2 3 4 


5 






t[2TT/U] 




1.0 
0.8 


- 




^/^X^^ 


- 


0.6 


- 




/ / 


: 


0.4 


- 




/ 


; 


0.2 


J 


'// 




; 


00 


J^ 


'y 




: 


C 


.00 


0.05 


0.10 0.15 0.20 


0.25 



VjU 



FIG. 1. (Color online) Top: Attenuation of first five revival os- 
cillations due to finite VjU = 0.005, 0.015, 0.025 (black, blue, 
red curves, V jU increasing from top to bottom), according to 
Q, !/ = 4, n = 2. Bottom: Attenuation factor of first coher- 
ence revival Av in ([5]) as a function of V jl] , for nv = 2, 4, 12 
(red, blue, black curves, nv increasing from right to left). 

The above result contains experimentally directly 
accessible information on the parameters U and V 
of the extended Bose-Hubbard model by measuring 
the quasimomentum distribution function rig(t) = 



^Y.^Ma, 



oiq^xi 



j^ ^_^i j\^i^ji^ - ~^i\ whose zero momentum q — Q 
part corresponds to the coherent fraction, where Xi are 
site centers. The height of the first few revivals at the in- 
stants 2'Kk/U , where k are positive integers, is highly sen- 
sitive to the value of V/U (assuming V <U), cf. Fig.[Tl 
where we plot the first order correlation function (setting 
Oi = 0), proportional to the zero momentum structure 
factor nQ{t), 



51 (i) = ^(«l"j) + ^■^■ 
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A significant attenuation of the revival peaks occurs al- 
ready for moderate values of V/U , which exceeds the 
experimentally observed damping in the V —^ Q (contact 
interaction) limit caused, e.g., by external trapping and 



finite temperature [^, [l^ . Hence the presently obtained 
homogeneous system result should readily be observable 
in experiments like [lOj, which aim at the minimization 
of external trapping (inhomogenity) effects, and can sen- 
sitively detect even small dipole strengths of particles 
trapped in an optical lattice. 

The attenuation of 2tt/U oscillations between revivals 
at 2'k/V increases with the product nv^ cf. Fig.[T] We 
take as a measure of the attenuation the value of Ay ~ 



I — {a[aj) /n at the first onsite revival tjj = 27r/[/, 



Ay = f — exp 
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The overdamped limit corresponds to Ay — >■ I. For small 
V/U, the attenuation is given by Ay = ni'{2TrV/U)^. 

We now compare our attenuation result ([5]) above 
with the real damping of the revival amplitude 
due to single-particle tunneling coupling between 
sites. In contrast to the small J result Aj = 
(2^)2[j,(2^_l)e-4nj2(2^) +j.e-3"i^|i^ij(n3)] (where 
IoiF{i,i} are modified Bessel and hypergeometric func- 
tions, respectively [la]), the small V result from (O is 
linear in i/ instead of quadratic. More importantly, the J 
damping decreases exponentially with filling n, while the 
V attenuation effect increases with the filling. However, 
while J/U will be necessarily small in the Mott phase (it 
is exponentially decreasing with increasing laser inten- 
sity), V/U can be of order unity [23|, so that the effect 
of V can be very pronounced while that of J will be sub- 
dominant, cf. Fig. [3] below. For U > V, the period of 
the first full-amplitude coherence revival is 2tt/V instead 
of 2tt/U. When U/V is not a rational number, no full 
height revival occurs, and so-called fractional revivals are 
obtained instead Q- 

Due to the long-range nature and slow fall-off of 
the dipolar interaction with distance ex 1/r'^, the next- 
nearest and next-to-next-nearest neighbors et cetera add 
other coupling terms of the form "^^ J2iij} ''^i^j to the 
Hainiltonian ([T]), where {ij} indicates a sum over next- 
nearest neighbors and A is their number, J2uk} = A at 
fixed k. Therefore, additional (slower) oscillation fac- 
tors will be added in ([3]). We take as a concrete exam- 
ple to illustrate this a 2D square lattice, with dipoles 
oriented perpendicular to the 2D plane. This geometry 
implies that the closest next-nearest neighbors are a dis- 
tance y/2a from a given site, such that the associated cou- 
pling Vnni/V = 1/V8, and the following next-to-next- 
nearest-neighbor sites at a distance 2a, giving a coupling 
Vnn2/V = 1/8. Therefore, while the revival at integer 
multiples of 27r/T/^ occurs with the full initial height when 
only nearest neighbors would be taken into account (i.e. 
for perfectly screened interactions at the NN distance), 
adding further next- nearest couplings, etc., makes the 
first full height revival peak occur much later. The full 
height revivals related to truly long-range interactions 
thus occur at a very late and hence not experimentally 
observable stage after the sweep. The influence of the 



Vnnh on the experimentally most relevant quantity, the 
attenuation of the first few revivals, is however small as 
long a.s V/U -^ 1. Our further discussion is based on the 
assumption that the Hamiltonian ([T]) with only nearest- 
neighbor couplings describes the system, keeping in mind 
the above limitation. 



III. MODIFYING THE GUTZWILLER 
APPROACH 

Conventionally, the equations of motion resulting from 
(fT|) are solved (numerically) with a Gutzwiller mean-field 
product state variational ansatz for the many-body wave- 
function, cf., e.g. [iiliii 
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with fi{ni, t ~ 0) ~ e '"• I' a"*/vni! for a coherent state 
© and where |ni, . . . ,71^/) = |ni) (g) • • • \ni) •••(g) {um) 
is the number basis state corresponding to noninteract- 
ing bosons occupying Wannier orbitals at each site. The 
above variational wavefunction ansatz is a bosonic ver- 
sion of the fermionic product state wavefunction origi- 
nally proposed by Gutzwiller [l|. 

The mean-field theory implicit in the Gutzwiller ansatz 
for the wave function is expected to be exact for the stan- 
dard onsite interaction Bose- Hubbard model when either 
n or ly [2l|, |30, |3l| are large; thus the "mean-field param- 
eter" nv, characterizing the validity of mean-field theory, 
should be sufficiently large. The Gutzwiller ansatz then 
becomes exact for contact interactions in infinite dimen- 
sion |32{ . Another extreme case of exact validity is a 
fully connected lattice [33|- The Gutzwiller wavefunc- 
tion, in the ground state and for contact interactions, 
has been argued to be qualitatively correct for lower, and 
hence easily experimentally realized, dimensions as well 
[3J]. The Gutzwiller ansatz ([6]) is however widely used 
for ground-state calculations also in the extended Bose- 
Hubbard model 0, |24| - [26J , where its rigorous validity for 
large mean-field parameter nu is much less obvious than 
in the onsite interaction case. As we will now argue, 
the extended Bose-Hubbard Hamiltonian cannot be de- 
coupled (in a nontrivial manner) into the sum of single- 
particle Hamiltonians, like the Hubbard model with con- 
tact interactions, where mean-field decoupling and the 
Gutzwiller ansatz have been shown to be equivalent [35| . 

The mean-field decoupling for the tunneling term, in 
the homogeneous case, = {at), leads from ([1]), when 
V^ = 0, to a sum of single-particle Hamiltonians, Hi = 

-Jiajcj) + (t)*a, - |0|2) + ^-hiih, - 1). On the other 
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FIG. 2. (Color online) Time evolution of gi with an initial 
coherent state for V/U = 0.05,0.1,0.3, from top to bottom, 
with filling n — 2. Solid line using the Gutzwiller variant 
((Tj), dotted blue line using the conventional Gutzwiller ansatz, 
Eq. p. 



hand, the nonlocal interaction term ex V, after apply- 
ing mean-field decoupling, just shifts the chemical poten- 
tial, and thus remains dynamically ineffective. In partic- 
ular, while the Gutzwiller wavefunction ([6]) reproduces 
the on-site-interaction modulations with period 27r/U, 
no modulations of gi with period 27r/I^ are created af- 
ter mean-field decoupling has been performed. Note that 
the (instantaneously) vanishing order parameter does not 
cause the failure of the conventional Gutzwiller wave- 
function ansatz. For local interactions, collapse and re- 
vival are captured accurately by the Gutzwiller wave- 
function, even though the mean-field amplitude vanishes 
when first-order coherence collapses. The inadequacy of 
the Gutzwiller ansatz for the time evolution of coher- 
ence in the extended Bose-Hubbard model thus is much 
stronger than that caused by the tunneling term ex J 
piSl \ll\ . The tunneling term can (at least in large di- 
mensions) be accurately decoupled in mean-field, while 
this is not possible in a nontrivial way for the V term in 
([1]), as argued above. We also observe in this regard that 
the expression ([3]) for the collapse and revival oscillations 
of _gi is nonperturbative in either 1/n or l/{ni>), so that 
there is no obvious mean-field limit of the exact quantum 
evolution of first-order coherence in deep optical lattices. 
The quantum evolution of coherence can therefore serve 
as a sensitive probe of various time-dependent ansatze 
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FIG. 3. (Color online) Time evolution for fixed V = 0.015 at 
two difTerent fillings n as indicated in the plots, and for vary- 
ing initial tunneling coupling J, J — (solid line, analytical 
result), J = 0.01 (red dashed) and J = 0.02 (blue dotted). 



for the many-body wavcfunction. 

Based on the above reasoning, we devise a variant of 
the dynamical Gutzwillcr wavcfunction, which takes into 
account that the time evolution according to ([1]) signif- 
icantly depends on the occupation of neighboring sites 
when V is appreciably large. We assume still a product 
of variational amplitudes fi in the many-body wavcfunc- 
tion, but arc now making the fi in addition dependent on 
the NN particle numbers with respect to the given site i. 



|GWV)=^...^if]/.(n,;,n[,,],i)l|ni,. 
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where nyi^ again indicates the set of all particle numbers 
on the NN sites j to a given i, cf. Eg.®. 

The ansatz ([7]) is physically corresponding to the fact 
that the on-site distribution function fi should also de- 
pend on the particle numbers in neighboring sites for 
finite V. This can be directly seen from the equations 
of motion for the variational functions following from 
the extended Bose- Hubbard model ([1]) when J = 0, see 
Eq. (jA.3|) in the appendix 
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dt 
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where /i is the chemical potential. Solving this equation 
generates the dynamical phase factors in ^. We have 
thus verified that as a benchmark for the validity of ([7]) 



the analytical result ([3|) is reproduced and hence the ex- 
act quantum evolution in the limit of zero tunneling. We 
illustrate this by comparing ([3]) with the evolution using 
either © or (O, see FiglH The evolution according to 
([7]) is identical with ([3]), while propagating the initial co- 
herent state with the conventional Gutzwiller ansatz (|6]) 
shows no modulation of the revival peaks with period 
27r/X^ at aU. 

We anticipate that the Gutzwiller variant (GWV) 
ansatz ([7]) will also be a good approximation of the dy- 
namical many-body state for a given t/, V when the tun- 
neling rate J is sufficiently small. We illustrate the effect 
of finite J by numerically computing in a ID lattice the 
time evolution of an initially coherent state, using the 
GWV tunneling energy functional (|A.7p derived in the 
appendix. We show the results in Fig. [31 There is the 
expected (real) damping occurring for finite J in addi- 
tion to the attenuation of the first few coherence revivals 
obtained by finite V . When nV > J (assuming U > V), 
which is the regime of validity for the ansatz ([7]) , the V 
attenuation will dominate the J damping; cf. Fig. [31 Con- 
versely, for J > nV, the conventional Gutzwiller ansatz 
(|6]) will again be applicable, i.e., yield in this regime a 
better approximation than the GWV ansatz ([7]). Fi- 
nally, coherent state amplitudes ^ apply in ((51) when 
one is entering deep into the superfiuid phase, which has 
J > nil, fiV. 



IV. NUMBER SQUEEZED INITIAL STATE 

To employ our variant of the Gutzwiller ansatz for a 
less simple initial state than the coherent state ^, we 
now investigate the NN coupling effect on collapse and 
revival when finite number squeezing is present in the 
initial state. We take as the initial state the numerically 
determined ground state of the Hamiltonian ([I} in one 
dimension for finite J /U and small system sizes [M = 7), 
by exact diagonalization (ED), so that the initial state 
is determined as the ground state of (jl} for a number- 
squeezed state. 

We evolve this number-squeezed initial state within 
the GWV wavefunction ansatz (O in a ID lat- 
tice. To implement GWV propagation, we pro- 
ceed by the following prescription. We con- 
vert the number basis amplitudes in the general 
many-body wavefunction obtained from ED, \^) = 
T,ni,...,nM /("-!' • ■ • ,"-A/)l"-i, ■ • ■,"-!/) by the prescription 
fi{ni) = 



Z-/ni ,....n.i 



-i:n.i-^-i,...,nM 



l/K,. 



,"m)P 



into real amplitudes of an initial product state of the 
GW type, and then propagate this initial state by the 
equations of motion of the GWV wavefunction, Eq. ([8]) . 
A comparison we made with the ED propagation of the 
full initial data shows that similar results arc obtained as 
long as ?7 = U /{vJ) does not become too large [i.e., larger 
than 0{n)\, in which case both initial coherenc es o i (0) 
are sufficiently close to the coherent state value [Sg . 
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FIG. 4. (Color online) Top: Time evolution of gi according to 
the GWV wavefunction ((7} for number-squeezed initial states 
with squeezing parameter 77 = Ui/{vJ) as indicated in the 
plot (z^ = 2, n = 3, |i — jl =3). The initial and final V — 0.03, 
final U = Uf = 1, as well as initial J — 1 and final J — are 
kept fixed, while rj is varied by changing the initial coupling 
U = Ui. Bottom: Time evolution for the GW ansatz ^, 
with identical parameter set and the same initial data. 



The basic result we obtain, upon increasing the squeez- 
ing parameter rj from zero (corresponding to a coherent 
state), is that revival peaks broaden but also that at the 
same time the attenuation of the first few [/-revivals is 
diminished, see Fig.|4] [33|- This is related to the fact 
that the narrower the (initial) distribution, the less phase 
spread in the dynamical phase factor is generated by 
the many-site coupling due to the V term in ([1]). In 
fact, the same broadening is obtained with the conven- 
tional Gutzwillcr ansatz (j6|), but like for the coherent 
initial state (77 = 0), no attenuation of the revival peak 
height is observed in distinction to the GWV evolution 
([7]), cf. Fig. m We note, in particular, that the Gutzwillcr 
approach also fails to reproduce that, for larger 77, the at- 
tenuation of the first revival peak heights relative to the 
initial coherence signal is reduced. 

Collapse and revival with reduced initial amplitude but 
less attenuation thus persists for not too large squeezing 
parameter 77. This opens up the possibility of probing 
the existence and time development of coherence revivals, 
and thus the interplay of quantum fluctuations and or- 
der parameter value, for various initial data sampled from 
around the transition point between supcrfluid and local- 
ized phases. The collapse and revival phenomenon ulti- 
mately owes its existence to the granularity of matter, 
i.e. the fact that particle number is quantized in integer 
multiples. Measuring the collapse and revival oscillations 
around the transition point with single-atom resolution 
|13 | should therefore afford interesting insights into the 



microscopic physics, that is the many-body wavefunction 
describing the transition. 



V. CONCLUSION 

The collapse and revival oscillations of first-order co- 
herence were shown to be a sensitive measure of the NN 
couplings in the extended Bose-Hubbard model. After 
a rapid quench from the superfluid side, the first few 
revival peaks are significantly attenuated for moderate 
NN couplings. Even relatively small dipole moments 
of molecules stored in optical lattices can thus be de- 
tected by the attenuation of the height of the first few 
revival peaks. The attenuation effect is strong enough to 
overcome masking inhomogenity effects like those caused 
by the external harmonic trapping, in distinction to the 
weak real damping caused by residual tunneling between 
deep wells in the localized phases. The Gutzwillcr variant 
ansatz we determined from the exact collapse and revival 
quantum evolution in the limit of zero tunneling is also 
potentially relevant for exploring the phase boundaries in 
the ground state of the extended Bose-Hubbard model. 
A particular example are supersolid phases in the limit 
of large NN coupling V, where a comparison should be 
instructive to those obtained by the standard Gutzwillcr 
approach |26l |. 

We expect the morphology of quantum phase revivals 
to be a valuable future tool for both the quantita- 
tive analysis of interaction couplings of Hamiltonians of 
atoms/molecules with long-range interactions stored in 
deep optical lattices, as well as a sensitive test-bed for the 
validity of various approximation schemes for the many- 
body wavefunction. 

Finally, a possible extension of the present work is to 
include pair-exchange terms between lattice sites, i.e., 
by adding a contribution Hpaii oc J2<ij>{^U^^ to the 
Hubbard Hamiltonian ([1]), see e.g. [3^, from which novel 
interaction-induced effects on the dynamical evolution of 
first- and second-order coherence can be expected |39| . 
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Appendix: Equations of motion from the GWV 
ansatz 



We delineate here the derivation of the equations of 
motion for the GWV wavefunction (O. For clarity and 
simplicity of notation, we restrict ourselves to a one- 
dimensional lattice; generalizations to higher lattice di- 
mensions are straightforward. 



We deduce the equations of motion variationally, start- 
ing from the functional 

£ = {GWY\^^t-Y,{H^-^^n,}\GWY), (A.l) 



where H^ == --'"I]j=j±i ^I'^j + f "j("-i - 1) + -j^i^t+i + 



2 ^i^i — 1 • 

We first consider the case with no tunneling, J — 0, 
for which ([7|) reproduces the exact quantum evolution. 
Due to the fact that the GWV wavefunction ^ is still 
a product state ansatz like the original Gutzwiller wave- 
function, the energy functional (jA.ip decomposes into 
contributions of a given site: 



£i= X! li{nii'n-i+ii'n'i-i,'t)fi-i{ni-i,'rii,ni-2,t)f*^i{ni+i,ni+2,ni,t) ■ ■ ■ 

X lidl''' - Hi{n^,ni+i,ni-i) + ^n,, /j(ni,nj+i, n,_i,i)/,_i(ni_i, ni,ri,j_2, t)/j+i(ni+i, ni+2, n,, t) ■ 



(A.2) 



where id^ now acts on fi{ni, rii+i, n^-i, t) only. Further- 
more, when J = 0, the phase factors of e.g. the functions 
fi+iini+i,ni+2,ni,t) and fi-i{ni-i,ni,ni-2,t), contain- 
ing all their time and NN particle number dependence, 
cancel with those of their conjugates in Si , so that they 
can be replaced by the onsite distributions of the initial 
state, fi±i{ni±i,ni,ni±2,t) -^ /°±i(»^i±i)- The equa- 
tions of motion then are 

[idi''-H,{ni,n,+i,n,^i,t) + fj,ni]fi{ni,ni+i,n,_i,t) = 0, 

(A-3) 
with the normalization condition 

J2ni \fii'n-i,ni+i,ni^i,t)\'^ = 1, leading to number con- 
servation in the form J2n "-jI/jI^-j: "-j+Ij ^i-i: i)P = "■■ 
The above equations of motion give the exact time 
evolution when J = 0. 

Due to [H{hi,ni±i), g{hi, ni±i)] = 0, where gini,hi±i) 
is an arbitrary operator function, a number of additional 
conserved quantities can be constructed. One such quan- 
tity is, choosing g{ni,n^±l) =hi+ fii+i + n^.i. 



(GWV| n, + n,+i + n,_i |GWV) 



(A.4) 



where the threefold product amplitude Fi is defined by 

Fi(ni,n,;+i,ni_i,i)) 

^^ 



[= /i(ni,ni+i,n,_i,t)/°_i(nj_i)/°+i(ni+i) for J = 
0]. Here, the Ji^i{ni±i^ni,t) are defined by 

stripping off the ni±2 dependence from the full 
/i±i amplitudes, e.g. /i_i(ni_i, ni,ni_2,i) = 

/i-i("-i-i, "-iii) exp[— i-jni_iri,i_2^]- This particular 
choice of g and the definition of Fi in the form above 
is motivated by C being conserved also in the J 7^ case 
(see below). In a lattice of any dimension and geometry, 
C = {v + l)n; for the ID lattice C = in. 



When J does not vanish and is small, the ansatz ([7]) 
can be applied to describe temporal evolution of the 
many-body state under certain conditions. Primarily, 
the latter consist in the requirement that the modulus 
of the variational amplitudes is not appreciably changed 
by finite and small tunneling and still dominated by the 
initial distribution, /f (n^), as well as that phase factors 
are to lowest order still given by the J = expression, 
exp [-^^{n^ni+l + n^n^-i)* - i^ni^ni - l)t] . 



The onsite contribution of the tunneling term £i_j 
to the functional (jA.l[) consists of four terms assigned 
with alai±i and ajj.^'^i- Consider one such term, 
— J (GWV| al_^_J^ai |GWV). Under the above assumption, 
to lowest order in J/V (assuming U > V), it can be 
shown to take the approximate form 



- J {al_^.iai) ~ -J ^ y/{ni + l)ni+iF*{Tii, ni+i,ni-i,t)Fi{ni + 1, Ui+i - 1, ni_i, t) ^ |/i+2(ni+2)pexp (iT/ni+2t) , 

(A.6) 



ni,ni±i 



ni+2 



I 

with similar expressions for the other three terms. In the at the end of the right-hand side can be omitted, so that 
short time domain t <^ lix jV we are interested in [cf. collecting all four terms, the tunneling contribution to 
the final paragraph of section II], the oscillating factor the functionals (|A.ip respectively (jA.2p becomes 



£i,.J - J ^ i^*(ni,nj+i,nj_i,t) \Jni{ni+i + l)Fj(nj - l^n^+i + l,ni_i,i) 



+ A/("i + l)'^i-l^j("i + l,7^^+l,7^^_l - l,t) 



where the normahzation condition now reads 
Sni.riiii l^i('^i'"-i±i'*)P = 1- Writing the interac- 
tion part as well in the product amplitudes (jA.Sp . 
the equations of motion can then be derived from the 
functional derivative d£i/dF*(ni,ni^i,ni^i,t) = 0. 



I 

Finally, one can easily verify that the tunneling term 
preserves the conservation of the sum of onsite and NN 
numbers, [-J{a]ai+i + ajai_i + al_^_^ai + aj_ja,),n£ + 
TT-i+i + hi^i] ~ 0, so that the conservation law (|A.4[) 
still holds, providing a benchmark for the accuracy of 
numerical calculations. 
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